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^ ■ Abstract 

The energies of the low- lying isoscalar and isovector udss configurations with spin-parity J p = 
+ , 1 + , and 2 + are calculated in a non-relativistic constituent quark model by use of the variational 
method. The contributions of various parts of the quark-quark interacting potentials including the 
s-channel interaction are investigated, and the effect of different forms of confinement potential is 
examined. The model parameters are determined by the same method as in our previous work, and 
they still can satisfactorily describe the nucleon-nucleon scattering phase shifts and the hyperon- 



3 
in 



nucleon cross sections. The parameters of the s-channel interaction are fixed by the masses of K 
and K* mesons, for which the size parameter is taken to be two possible values. When it is chosen 
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\ as the same as baryons', the numerical results show that the masses of all the udss configurations 



are higher than the corresponding meson-meson thresholds. But when the size parameter for the 
K and K* mesons is adjusted to be smaller than that for the baryons, the udss configuration with 
1 = and J p = 1 + is found to lie lower than the K*K* threshold, furthermore, this state has a 
very small KK* component and the interaction matrix elements between this state and KK* is 
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comparatively small, thus its coupling to the KK* channel will consequently be weak and it might 
be regarded as a possible tetraquark candidate. 
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I. INTRODUCTION 



So far, all the observed hadrons can be classified into two types, i.e. "baryons" composed 
of qqq and "mesons" composed of qq. But in principle, the QCD fundamental theory doesn't 
exclude the existence of the states containing more than three quarks, i.e. the so-called 
multi-quark states. Since Jaffe predicted the H particle (uuddss) in 1977 the research 
on multi-quark states has always been an attractive topic for nearly three decades in both 
theoretical and experimental studies. But up to now, there is no convincing evidence of their 
existence in experiments. In 2003, LEPS Collaboration reported the possible existence of 
ihe + pentaquark and subsequently some laboratories also reported the similar results 
3]]. At the same time, several laboratories reported the negative results 3J. Although its 
existence is still questioned, the + particle has motivated a number of theoretical and 
experimental studies of pentaquarks and further the multi-quark states. 

Besides dibaryon and pentaquark, the possible udss tetraquark is another interesting 
multi-quark system, and many works have been devoted to the investigation of this state 
in the past few years [l, O, [g, Q, [s, [9], Q. In Ref. 4], using the MIT bag model, Jaffe 
performed a wide study for the spectrum of the Aq states, and he predicted the masses of 
the isovector J p = + and isoscalar J p = 1 + udss states to be 1.55 GeV and 1.65 GeV, 
respectively. I, Ref. & V ijaIld e et al. analyzed the ^ system for both is o 8 p ln / = 
and 1 = 1 channels in a constituent- quark model, and they didn't find any stable udss 
tetraquark state. In Ref. [6j, Burns et al. claimed that the @ + particle suggest there 
should exist a udss tetraquark state with J p = l - and mass around 1.6 GeV. This state has 
a strong color- magnetic attraction and decays into KK channel via P-wave with a width 
around 10 — 100 MeV. In Ref. |7j], Karliner and Lipkin argued that the isoscalar J p = + 
udss tetraquark state is a cousin of the + pentaquark, and for this state the lowest allowed 
decay mode is a four-body KKwk channel with a very small phase space and a distinctive 
experimental signature. However, we should note that from the non-clustered quark degree 
of freedom, the isoscalar J p = + udss state within the spatially symmetric configuration 
is not allowed due to the Pauli principle. In Ref. [8], Kanada-En'yo et al. studied the udss 
system in the framework of the flux-tube quark model, and they pointed that the isoscalar 
J p = 1 + udss is a stable and low- lying tetraquark state with a mass around 1.4 GeV. 
It can decay into K*K via S-wave with its width around 20-80 MeV, while the isovector 
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J p = + udss tetraquark does not exist. Recently, Cui et al. calculated the masses of 
the isoscalar J p = 1 + udss systems by using the color- magnetic interaction Hamiltonian 
with SU(3) flavor symmetry breaking, and they found the udss tetraquark lies around 1347 
MeV with a narrow width since it can not decay into the KK* channel. In Ref. [k]], Chen 
et al. performed a QCD sum rule study for the isovector J p = + udss system, and they 
obtained a tetraquark state with a mass around 1.5 GeV. 

To sum up, the existence and properties of the possible udss tetraquark state are presently 
model dependent. Further theoretical and experimental investigations of this state via dif- 
ferent approaches seem to be significant and essential. 

It is a general consensus that QCD is the underlying theory of the strong interaction. 
However, as the non-perturbative QCD effect is very important for light quark systems 
in the low energy region and it is difficult to be seriously solved, people still need Q CD- 
inspired models to be a bridge connecting the QCD fundamental theory and the experimental 
observations. In the past few years, we have developed a non-relativistic constituent quark 
model, which has been quite successful in reproducing the energies of the baryon ground 
states, the binding energy of the deuteron, the nucleon-nucleon (NN) scattering phase 
shifts, and the hyperon-nucleon (YN) cross sections [llj. In this model, the quark-quark 
interaction contains confinement, one-gluon exchange (OGE) and boson exchanges stemming 
from scalar and pseudoscalar nonets, with the boson exchange potentials deduced from a 
linear interacting Lagrangian which is — under the ch,a, —.nation Q, and 
the short range quark-quark interaction is provided by OGE and quark exchange effects. 

Actually it is still a controversial problem in the low-energy hadron physics whether 
gluon or Goldstone boson is the proper effective degree of freedom besides the constituent 
quark. Glozman and Riska proposed that the Goldstone boson is the only other proper 



effective degree of freedom [13j, ll4j. But Isgur gave a critique o 



the boson exchange model 



161 ]. Anyway, it is still an 



and insisted that the OGE governs the baryon structure [15, 
open problem in the low-energy hadron physics whether OGE or vector-meson exchange 
is the right mechanism for describing; the short-range quark-quark interaction, or both of 



them are important. Thus in Ref. 



171 ] we further extended our original constituent quark 



model to include the coupling of the quark and vector meson fields. The OGE that plays 
an important role in the short-range quark-quark interaction in our original model is now 
nearly replaced by the vector boson exchanges. This model has also been successful in 
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reproducing the energies of the baryon ground states, the binding energy of the deuteron, 
and the nucleon-nucleon scattering phase shifts 



Recently, we have extended our constituent quark model from the study of baryon 



scattering 
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processes to the baryon-meson systems and the pentaquark + state 
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baryon 
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241 ]. We found that some results are similar to those given by the chiral 



unitary approach study, such as that both the AK system with isospin 1=1 and the UK 



system with / = 1/2 have quite strong attractions 
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scattering 
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2l|. In the study of the KN 



19] , we get a considerable improvement not only on the signs but also on 



the magnitudes of the theoretical phase shifts comparing with other's previous work. We 



also studied the structures of the + particle 23j, |24J , and found that the theoretical mass 



of + is much higher than the experimental value, thus we concluded that either the G + 
particle does not exist or it can not be explained by a five-quark cluster, which is consistent 
with the high-statistic experimental negative results. 

In this paper, using the variational method, we study the structures of the isoscalar and 
isovector udss configurations with spin-parity J p = + , 1 + , and 2 + in our non-relativistic 
constituent quark model with meson and one gluon exchanges included. The contributions 
of various parts of the quark-quark interacting potentials including the s-channel interaction 
are investigated. The effect of different forms of confinement potential is examined, and the 
configuration mixing between the states with same quantum numbers is considered. Th e 
model parameters are determined by the same method as in our previous work {20 . 



21 



22], 



and they still can satisfactorily describe the nucleon-nucleon scattering phase shifts and the 
hyperon-nucleon cross sections. The s-channel quark- ant iquark interaction is complicated, 
and we take the zero momentum approximation and fix the corresponding parameters by 
fitting the masses of K and K* mesons. With the size parameter of K and K* taken to be the 
same as baryons', the numerical results show that the masses of all the udss configurations 
are higher than the corresponding meson-meson thresholds. But when the size parameter 
for the K and K* mesons is adjusted to be smaller than that for the baryons, the udss 
configuration with 1 = and J p = 1 + is found to lie lower than the K*K* threshold, 
furthermore, this state has a very small KK* component and the interaction matrix elements 
between this state and KK* is comparatively small, thus its coupling to the KK* channel 
will consequently be weak and it might be regarded as a possible tetraquark candidate. 
The paper is organized as follows. In the next section the framework of the non-relativistic 
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constituent quark model we used and the wave functions of the udss configurations are briefly 
introduced. The calculated energies of six udss states are shown in Sec. Ill, where some 
discussion is presented as well. Finally, the summary is given in Sec. IV. 



II. FORMULATION 



A. Model 



Our non-relativistic constituent quark model has been widely described in the literature 
and the details can be found in these references. Here we just give the salient 
features of our model. 

The total Hamiltonian of the udss systems in the model can be written as 

h = J2t 1 -t g + v 12 + v m +J2 v& (!) 



i=l,2 
.7=3,4 



where Tq is the kinetic energy operator for the center-of-mass motion, and Vu, V34 and Vq 
represent the quark-quark, antiquark-antiquark and quark-antiquark interactions, respec- 
tively, 

Vl2 = \/° GE + + (2) 
where V^ GE is the OGE interaction, 



12) 



1 14 1 

— 2 + — 2 + Q (°1 ' °2 

mf m 2 6 m,\m2 



(3) 



A— 1 " [r 12 2 

and the confinement potential ? , instead of the quadratic form as used in our previous 
work, is taken as the linear one, 

= - (A? • A c 2 ) (a c 12 r 12 + <) . (4) 

represents the effective quark-quark potential induced by the one-boson exchanges. 
In our original constituent quark model, Vi 2 includes the scalar boson exchanges and the 
pseudoscalar boson exchanges, 

8 8 
a=0 a=0 
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and when the model is extended to include the vector boson exchanges, can be written 

as 

8 8 8 

v i2 = E ^.( p «) + E ^.( p «) + E v p*to)- ( 6 ) 

a=0 a=0 a=0 

Here <t , a 8 are the scalar nonet fields, 7r , .., 7r 8 the pseudoscalar nonet fields, and p , .., p 8 
ihe vector nonet fields. The expressions of these potentials can be found in the literature 



]J, 1181,1191,0 • 



V34 in Eq. (1) represents the antiquark-antiquark interaction, 

Vu = Vi GE + Vit + VH (7) 

where Vg GE and V^ onf can be obtained by replacing the A^ ■ A£ in Eqs. (3) and (4) with 
A|* • Af*, and V^ 1 has the same form as V$. 

Vfj in Eq. (1) represents the quark-antiquark interaction, 

Vq = Vf GE + V t f nf + V$ + V%, (8) 

where V| OGE and V^ onf can be obtained by replacing the A^ ■ A| in Eqs. (3) and (4) with 
— A^ • AS*, and V^ h can be obtained from the G parity transformation: 

vi = E(- i ) Gfc ^ h ' fc > ( 9 ) 

k 

with (— l) Gk being the G parity of the kth meson. V% denotes the s-channel quark-antiquark 
interaction. For the udss system, the s-channel interaction includes K and K* exchanges, 

v$=v*+vr- (io) 

Taking the zero momentum approximation, the spatial part of Vf. can be expressed as a 
delta-function. To flatten the delta-function, we replace it by the Yukawa function, then 
V S K and V S K * can be expressed as: 

v? = c -< (^p)_ (i±M)^ (2)^-, („, 

and 

vr=c r (31^^(1+^^)^ (12) 

where C K and C K * are treated as parameters and we adjust them to fit the masses of K 
and K* mesons. 
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B. Parameters 



TABLE I: Model parameters. The meson masses and the cutoff masses: m a / = 980 MeV, m K = 980 
MeV, m e = 980 MeV, = 138 MeV, m K = 495 MeV, m v = 549 MeV, m. q , = 957 MeV, m p = 770 
MeV, m K , = 892 MeV, = 782 MeV, ms = 1020 MeV, and A = 1100 MeV. 





Model I 


Model II 

j ciiv / tf cnv 


Model III 

fnhir/Clnhir = 2/3 
J CIIV / iJ CIIV j 


\ J " l " lJ -/ 


0.5 


0.45 


0.45 


m u (MeV) 


313 


313 


313 


m s (MeV) 


470 


470 


470 


q 2 


0.766 


0.056 


0.132 


q 2 

Us 


0.846 


0.203 


0.250 


9ch 


2.621 


2.621 


2.621 


9chv 




2.351 


1.973 


m a (MeV) 


595 


535 


547 


a c uu (MeV/fm) 


87.5 


75.3 


66.2 


a c us (MeV/fm) 


100.8 


123.0 


106.9 


a c ss (MeV/fm) 


152.2 


226.0 


196.7 


a* (MeV) 


-77.4 


-99.3 


-86.6 


at (MeV) 


-72.9 


-127.9 


-109.6 


af s (MeV) 


-83.3 


-174.20 


-148.7 



The harmonic-oscillator width parameter b u is taken to be 0.50 fm in our original con- 
stituent quark model, and when the vector boson exchanges are included, b u is taken to 
be 0.45 fm. This means that the bare radius of baryon becomes smaller when more meson 
clouds are included in the model, which sounds reasonable in the sense of the physical pic- 
ture. The up (down) quark mass m u ^) and the strange quark mass m s are taken to be the 
usual values: m u (d) =313 MeV and m s = 470 MeV. The coupling constant for scalar and 
pseudoscalar meson field coupling, g ch , is determined according to the relation 

9ch _ (3\ 9NNn m u /-|o\ 

4tt ~ \5J 4tt M 2 n ' 1 ' 

7 



with the empirical value g 2 NNlv /Air = 13.67. The coupling constant for vector coupling, of the 



vector- meson field is taken to be g c h v = 2.351, the same as used in the NN case [17j|. The 
masses of the mesons are taken to be the experimental values, except for a meson. The mo- 
is adjusted to fit the binding energy of the deuteron. The OGE coupling constants and the 
strengths of the confinement potential are determined by baryon masses and their stability 
conditions. All the parameters are tabulated in Table I, where the first set is for our original 
constituent quark model, the second and third sets are for the models with vector meson 
exchanges included by taking f c hv/9chv as and 2/3, respectively. Here / c h v is the coupling 
constant for tensor coupling of the vector meson fields. 

From Table I one can see that in models II and III, g\ and g 2 are much smaller than 
the values in model I. This means that when the coupling of quarks and vector meson fields 
is included in the non-relativistic constituent quark model, the coupling constants of OGE 
will be greatly reduced. Thus the OGE that plays an important role of the quark-quark 
short-range interaction in our original constituent quark model is now nearly replaced by 
the vector-meson exchange. In other words, the mechanisms of the quark-quark short-range 
interactions in these models are quite different. 

We'd like to mention that our previous work concentrated on the hadron-hadron inter- 
actions and it can be strictly proved that different forms of confinement potential does not 
make any visible influence on the theoretical results since the two hadrons are treated as two 
color-singlet clusters. In this work we adopt a color linear confinement potential to study 
the udss one-cluster system, and the method of parameters determination is the same as in 
our previous work. Naturally, the three sets of parameters in Table I still can satisfactorily 
describe the energies of the baryon ground states, the binding energy of the deuteron, the 
nucleon-nucleon scattering phase shifts and the hyperon-nucleon cross sections. 
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C. udss configurations 



According to the Pauli principle in each quark (antiquark) pair, in the spatial symmetry 
case only six configurations are permitted for the udss system: 

|1) = ({ud}l{ss}l) 



o,j p = r 



|2> = ({«d}g{55}«) 

|3> = (Ho{s«}i)i 
|4) = (M]?{^) 1 

I=l,J P = l + =► \5) = ({ud}l{ss}l) 1 



J=1,J P = 2+ =► | 6 ) = ({ud}l{ss}l) 2 

The s quark has isospin-zero so the total isospin of the udss configuration state is deter- 
mined by the isospin of the u, d quarks. In the above expressions, { } and [ ] represent the 
flavor symmetry and antisymmetry, respectively, and the superscript is the representation 
of the color SU(3) group, the subscript is the spin quantum number. Making a re-coupling 
calculation, we can express the above six states into two quark-antiquark pairs, including 
two color octet qq pairs and two color singlet qq pairs with KK, KK* and K*K* quantum 
numbers. The corresponding expressions are given as follows: 

|1> EE ({«*}?{55}?) o = \ ((us)l(ds)l) l - yj((«5)i(d5)i); 

-\J\ ((usTMs)*)] + yjl ((us)l(ds)l)l, (14) 

|2> ee ({«d}g{w}j) o = ((us)l(ds)l)l + yi ((us)\(ds)l)l 

+\[^ (Mg(ds)g); + i ((«5)?(d5)?); , (15) 

is) ee (hJ{5s}?) i = yj((«5)5(d5)a; - + ((«5)i(dg)a; 

+ yf (Ms^s; - y| ((«s)?(d5)S); , (i6) 

|4) EE (H?{55}S) i = -y|((«5)5(d5)i); + y|((«5)i(d5)5); + y|((«5)i(df)i); 



|5> ee [{ud}l{ss}l) i = ^{(us) 1 (ds)\)\ + ^Jl((us){(ds)l)l 

-y|((«5)g(d5)S); -^Jl((us)l(ds) s )l (18) 

|6) ee ({«d}?{55}?) 2 = ^1 ((«5)i(d5)i);- (19) 

where (us)(ds) represents 072 [(us)(ds) + (ds)(«s))] for Eqs. (14), (15), (18) and (19), 
and denotes ^l/2[(us)(ds) - (ds)(us)] for Eqs. (16) and (17). 

In the actual calculation, the configuration mixing between the states |1) and |2), as well 
as |3) and |4) is considered since these states have same quantum numbers. 

III. RESULTS AND DISCUSSIONS 



TABLE II: Energies (in MeV) of udss six configurations in our constituent quark models (without 
configuration mixing). 



Configurations 


Model I 


Model II 


Model III 


Threshold 


J p = 0+ 










({ud}l{ss}l) 


1679 


1650 


1650 


KK (990) 


({«d}g{is}jj) o 


1828 


1865 


1833 


KK (990) 


J p = 1+ 










(Mg{5s}?) i 


1698 


1704 


1698 


KK* (1387) 


(Mf{o}») i 


1803 


1812 


1798 


KK* (1387) 


(W}f{,- s -}?) i 


1765 


1750 


1742 


KK* (1387) 


J p = 2+ 










({nd}?{ S - S "}?) 2 


1904 


1906 


1885 


K*K* (1784) 



We calculate the energies for six low configurations of udss system in our non-relativistic 
constituent quark models. The model parameters we used are shown in Table I, which can 
reproduce the NN phase shifts reasonably. As mentioned above, the parameters of the s- 
channel interactions are fixed by fitting the masses of K and K*. For reducing the input 
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parameters, here the size parameter for K and K* is taken to be the same as that for 
baryons. The calculated results (without configuration mixing) are given in Table II. From 
this Table, we can see that the results from model I, II and III are quite similar, although 
the mechanisms of the quark-quark short-range interactions are different in these models, 
i.e., one is from OGE and the two others are from vector meson exchanges. This means that 
the OGE and the vector- meson exchange can give similar contributions in the udss system, 



the same situation as in the nucleon-nucleon and kaon-nucleon systems 1?], 

In Ref. {9], using the color-magnetic interaction Hamiltonian with SU(3) flavor symmetry 
breaking, Cui et al. argued that the strong attractive color-magnetic interaction can reduce 
the energies of the udss systems, and they found a / = and J p = 1 + udss tetraquark state 
with a mass around 1347 MeV. In our constituent quark model calculations, in model I, the 
color-magnetic interactions are attractive in both the isovector J p = + ({ud}i{ss}fj state 
and the isoscalar J p = 1 + ([wi]o{ss}f) 1 state, while they are repulsive in the other four 
configurations, and in models II and III, the OGE is largely reduced and the color-magnetic 
attractions are almost replaced by p exchange. Furthermore, in all the models, the a and 
Ti exchanges provide more attractive interactions in the isovector J p = + ({ud}l{ss}fj Q 
state and the isoscalar J p = 1 + ([wi]o{ss}f) 1 state than in the other four configurations. 
Thus the energies of the [{ud}l {ss}f) and ([wc?]o{ss}f) x states are respectively the lowest 
one in J p = + and J p = 1 + cases in various models. However, due to the high kinetic 
energies, the attractive interactions are not strong enough to reduce the energies of these 
two states to be lower than the corresponding meson-meson thresholds (see Table II). 

The states ({ud}\{ss}l) Q and ({ud} 6 {ss}l) Q in J p = 0+ well as ([Mrf]o{ss}i) 1 

and ([mc?]5{ss}q) in J p = 1 + case have the same quantum numbers. Thus the configuration 
mixing between them has to be considered. The results are shown in Table III. Comparing it 
with Table II, one can see that the configuration mixing effect is significant, it can shift the 
energies over 80 MeV for most configurations. In the isovector J p = + case, now the energy 
of the lowest state is reduced about 170—290 MeV in various models, and in the isoscalar 
J p = 1 + case, it is about 80—130 MeV. Even though, the energies of these configurations 
are still higher than their corresponding meson-meson thresholds. In other words, the stable 
udss tetraquark state cannot yet be obtained. 

The confinement potential is phenomenological, and usually it is taken as linear, quadratic 
or error function form. Here we consider these three various forms of the confinement 
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TABLE III: Energies (in MeV) of udss states with configuration mixing considered. 



Model I Model II Model III 

I = 1, J p = 0+ 1394 1483 1482 

1963 1952 1924 

I = 0,J P = 1+ 1567 1619 1615 

1887 1878 1861 

I = 1, J p = 1+ 1765 1750 1742 

I = 1, J p = 2+ 1904 1906 1885 



potential to see the corresponding effects. The results, from our original non-relativistic 
constituent quark model, are shown in Table IV. In this Table, r, r 2 and erf represent the 
confinement potential adopted as linear, quadratic and error function form, respectively. We 
can see that the energies with the confinement potential taken to be the error function form 
are always the lowest. But the difference of various confinement potentials is less than about 
30 MeV, which means different form of the confinement potential has no significant effect 
on the energy of the udss state. 

TABLE IV: Energies (in MeV) of the udss states in our original constituent quark model with con- 
figuration mixing and three different forms (linear, quadratic and error function) of the confinement 



potential considered. 

r r 2 erf 

/ = 1, J p = 0+ 1394 1402 1390 

1963 1996 1945 

/ = 0, J p = 1+ 1567 1568 1567 

1887 1907 1876 

/ = 1, J p = 1+ 1765 1770 1763 

/ = 1, J p = 2+ 1904 1920 1895 



As we know, the s-channel quark- ant iquark interaction mechanism is a complicated and 
unclear problem. In the study of the structure of the + particle |23j, it has been pointed 
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TABLE V: Energies (in MeV) of udss mixing states without and with s-channel interaction in 
different models. 



without s-channel interaction with s-channel interaction 





I 


II 


III 


I 


II 


III 


I = 1,J P = 0+ 


1699 


1910 


1890 


1394 


1483 


1482 




2041 


2001 


1991 


1963 


1952 


1924 


I = 0, J p = 1+ 


1749 


1825 


1828 


1567 


1619 


1615 




1984 


1964 


1960 


1887 


1878 


1861 


I = 1, JP = 1+ 


1876 


1922 


1912 


1765 


1750 


1742 


/ = 1, JP = 2+ 


1964 


1939 


1935 


1904 


1906 


1885 



out that how to treat the s-channel interaction reasonably is very important. In the udss 
system, the effect of the s-channel interaction should also be examined. We completely omit 
this interaction to see its influence, and the results are shown in Table V, compared with 
those with s-channel interaction. One can see that the s-channel interactions offer quite 
strong attractions in all these six udss configurations and thus can reduce the energies of 
these states about several hundreds MeV in all the models. In this sense, the effect of the 
s-channel interactions is significant and un-negligible in the udss system. 

TABLE VI: Energies (in MeV) of udss states with the configuration mixing considered and the 



size parameter of K, K* taken as 0.4 fm. 

Model I Model II Model III 

I = 1, J p = 0+ 1602 1573 1572 

1857 1909 1882 

J = 0, J p = 1+ 1577 1623 1618 

1768 1833 1817 

/ = 1, J p = 1+ 1771 1754 1745 

I=l,J p = 2+ 1821 1872 1852 
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In the above calculations, the size parameter of K and K* mesons are taken to be the 
same as the baryons' in order to reduce the free parameters. It seems more reasonable to 
choose the size parameter for the mesons to be smaller than that for the baryons. We then 
perform a variational method calculation for the energies of the udss states with the size 
parameter of K and K* taken to be 0.4 fm, and the results are shown in Table VI. We 
notice that one of the / = J p = 1 + udss configurations is a very interesting state in model 
I. Its energy is 1768 MeV, lower than the threshold of K*K*, and the corresponding root 
mean square radius of this state is about 0.57 fm. Especially, the structure of this state 
is very interesting since it contains very few components of KK*, as can be seen from the 
expression of its wave function given below: 

|4>'= - 0.17 ((us)l(ds)l)l + 0.17 ((uS)\(d8)l)\ + 0.71 {(uS)\(ds)\)\ 

- 0.47((«s)g(ds)?)J + 0.47 ((us)l(ds)*)\ + 0.0056 ((us)?(tfc)?)J . (20) 

From this equation it is clear to see that besides 43.8% part of two color octet qq pairs, the 
component of two color singlet qq pairs is 50.4% for K*K* and only 5.8% for KK*, which 
means the K*K* component is dominate and comparatively the KK* component is very 
small in the color-singlet qq-qq part, thus the / = J p = 1 + udss configuration might have 
a few possibility decaying into K and K*. 

Furthermore, in order to see the effect of the coupling between this / = J p = 1 + udss 
state and the KK* channel, we calculate the interaction matrix elements between the |4)' 
state and KK* channel in model I by treating K and K* as two clusters with a distance S, 

(V) = (4'\J2^\KK*^(R-S)), (21) 

j£K* 

— * — * 

where ip(R — S) is the relative motion wave function of the two clusters K and K*, and for 
simplicity we take it as 

ip(R -S) = (ccvVvr) 374 exp[-ua{R - S) 2 /2\, (22) 

here S is the generator coordinate which can qualitatively describe the distance between the 
two clusters. The calculated results are shown in Figure 1, along with the contributions of 
((us)l(ds)\)\, ((us)\(ds)\)\ and the two color-octet qq pairs components in the |4)' state. 
In Fig. 1, M denotes the interaction matrix elements, and S can qualitatively describe the 
distance between the two clusters K and K*. 
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FIG. 1: The matrix elements of the interaction potential. The solid line represent the matrix 
elements of interaction potential between the |4)' state and the KK* channel. The dotted, dash- 
dotted and dashed lines represent the contributions of ((us)g((is)})*, (^us)\{ds)\) v and the two 
color-octet components in 4)' state, respectively. 

From Fig. 1 one can see that the contribution of color-octet components to the interaction 
matrix elements is the smallest which is apparent and can be easily understood, while the 
((us)\(ds)Y) x component also has a small contribution due to the small component in |4)' 
(smaller than 6%), although the matrix element {{{us)\{ds)\)\ \ <eK Vij\KK*ip(R — S)) 
is comparatively big. The solid line in Fig. 1 clearly tells us that the total contribution of 
4)' to the interaction matrix elements, i.e. (4'| Vij\KK*ip(R — S)), is very small, less 

than 6 MeV and inclines to zero with the increase in distance between K and K*. Hence, the 
coupling between this interesting |4)' state and the KK* channel can be regarded as quite 
small. This means the / = J p = 1 + udss state has a few possibility decaying into two 
separate K and K* via interaction potential. Furthermore, since the energy of this state, 
1768 MeV, is lower than the K*K* threshold (1784 MeV), it cannot decay into K*K* final 
state. This means this state would possibly have a narrow width, and might be treated as 
a good candidate for the udss tetraquark state. 

As discussed in Ref. [9|, it seems worth searching this state in K + d — » p + p + K~ + T + 
or J/i/)(Y) — > K~ + K° + T + channels in the future experiments [Here T + denotes the 
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udss tetraquark state with 1 = and J p = 1 + ]. The experimental information about the 
existence of this state will help us test the validity of the application of our constituent 
quark models to the study of multi-quark states, although these models are quite successful 



in the investigations of NN, NY, and KN interactions 



IV. SUMMARY 

The structures of udss states with J p = + , 1 + , and 2 + are studied in our non-relativistic 
constituent quark models with the meson and one gluon exchanges included in the quark- 
quark interaction potentials. We calculate the energies of six low-lying udss configurations 
by use of the variational method. The configuration mixing between the states with same 
quantum numbers are considered. The effect of different forms of color confinement poten- 
tials and the contributions of s-channel qq interactions are also examined. The results show 
that the different forms (linear, quadratic, and error function) of the confinement potential 
just give similar contributions, and the s-channel interactions can reduce the energy of the 
udss system several hundred MeV. With the model parameters determined by the same 
method as in our previous work, the calculated energies of all the udss configurations are 
higher than the corresponding meson-meson thresholds. But when the size parameter for 
the mesons is adjusted to be 0.4 fm, a value smaller than that for the baryons, the udss 
configuration with 1 = and J p = 1 + is found to lie lower than the K*K* threshold, fur- 
thermore, this state has a very small KK* component and the interaction matrix elements 
between this state and KK* is comparatively small, thus its coupling to the KK* channel 
will consequently be weak and it might be regarded as a possible tetraquark candidate. A 
further dynamical calculation would be done in the future work. 
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